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ABSTRACT 

The astrophysical p process, which is responsible for the origin of the proton rich stable nuclei heavier 
than iron, was investigated using a full nuclear reaction network for a type II supernova explosion 
when the shock front passes through the O/Ne layer. Calculations were performed with a multi-layer 
model adopting the seed of a pre-explosion evolution of a 25 solar mass star. The reaction flux was 
calculated to determine the main reaction path and branching points responsible for synthesizing the 
proton rich nuclei. In order to investigate the impact of nuclear reaction rates on the predicted p- 
process abundances, extensive simulations with different sets of collectively and individually modified 
neutron-, proton-, a-capture and photodisintegration rates have been performed. These results are 
not only relevant to explore the nuclear physics related uncertainties in p-process calculations but are 
also important for identifying the strategy and planning of future experiments. 
Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: general 



1. INTRODUCTION 

The majority of the observed heavy nuclei above iron 
have been produced by neutron induced nucleosynthesis 
processes such as the slow neutron capture process (s 
process) in low mass AGB stars and massive red giant 
stars and the rapid neutron capture process (r proc ess) 
in the supernova shock front ijKappeler et al. |l!989|) . 

In nature, 35 nuclei can be found on the neutron de- 
ficient side of the valley of stability ranging from 74 Se 
to 196 Hg, which are shielded against production by neu- 
tron capture processes. For the nucleosynthesis of these 
isotopes the p process has been proposed as the most 
likely scenario. The nature of the p process is still un- 
der debate and it is possible that it represents differ- 
ent independent nucleosynthesis scenarios. The contri- 
bution from charged particle induced nuclear reactions 
are most likely negligible because the high Coulomb bar- 
rier reduces the associated reaction rates significantly, 
though the possibility of some contributio n to the light- 
est p - nuclei cannot be entire ly excluded l|Schatz et al. I 
119981 Ijordan fc Mever Il2004|) . In its current interpreta- 
tion the p process is described as 7 induced photodis- 
integration of stable nuclei in the shock front of type 
II su pernovae l|Wooslev fc Howard I Il97l iRavet et all 
1990) or alternatively in the deflagration flam e of a type 
I supernova detonation l|Howard et al. 1119911) . More re- 
cently pre-explosive sites for the p-process have been 
suggested in the O-Ne burning zone of massive stars. 
This feature emerges in one dimen sional stellar evolu- 
tion models l)Rauscher et al. ll2002ft and also in two di- 
mensional simulations of convective oxygen rich burn- 



ing zones. However, in the two dimensional study the 
calculation of the associated nucleosynthesis pattern has 
been limited to a one dimensional model adaptation and 
does not yet include quan titative abundance predictions 
ijBaleisis fc Arnett1l200H) . Supercritical accretion disks 
associated with jets in supe rnovae have also been pro- 
posed as a possible scenario l)Fuiimoto et al. 1120031) . 

The bulk of p isotopes represents a small fraction of the 
total abundance and is presently bas ed on the analysis of 
meteorite data l) Anders et al. l fl989). This is illustrated 
in Fig. ^ where the abundance distributions are com- 
pared for all those heavy isotopes, which can be entirely 
ascribed to a specific production process. In most cases 
the predicted p-process abundances agr ee within a fac- 
tor of three with the observed v alues l)Lambert 1119921 
iMever Ill99l IRavet et al. Ill995|) . However, there still 
are significant discrepancies for the light p nuclei 92,94 Mo 
and 96,98 Ru, which are largely overabundant compared 
to the model predictions, as well as for the Dy and Gd 
p- isotopes in mass region A = 150 — 160. Such discrep- 
ancies are reviewed in more detail in lArnould fc Gorielv I 
(2003), who give a comprehensive review of the present 
interpretation an d questions related to p-pr ocess nucle- 
osynthesis. While[Arnould & Gorielv (2003) discuss the 
global aspects of p-process nucleosynthesis, the present 
work concentrates on the nuclear physics related uncer- 
tainties in p-process model predictions and in particular 
on the identification of critical p-process reaction rates, 
which are most important for further experimental stud- 
ies. 

The presently favored scenario for the p process are 
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type II supernova (SN) expl osions ijArnould fe GorieTvl 
20031 lHavakawa et al. 112001 . The emerging shock front 
causes a rapid increase in temperature and density in the 
different layers of the pre-supernova star. As displayed 
in Fig. |3 peak temperatures between 1.7<Tg <3.3 are 
reached in the Ne/O layer of the presupernova envelope 
ijRavet et al. Ill995j) . The associated intense photon flux 
induces a range of photo-disintegration processes shifting 
the existing distribution of seed abundances to the pro- 
ton rich side of the valley of stability by (7,n)-reactions. 
When this process becomes less efficient because the neu- 
tron binding energy increases with neutron deficiency, 
the reaction flux is maintained by the (7,p) and (7, a) 
channels. After the shock front has passed the Ne/O 
layers, temperature and density drop exponentially and 
the unstable proton rich nuclei decay back to the valley 
of stability. In chapter 3 the p-process reaction flux will 
be discussed more quantitatively in terms of our present 
model predictions. 

The description of the entire synthesis process for the p 
nuclei requires a comprehensive reaction network involv- 
ing far more than ten thousand reactions. With very few 
exceptions the astrophysical reaction rates have been cal- 
culated by means of the statistical Hauser-Feshbach (HF) 
model, which can only be applied for nuclei with high 
level densities. The statistical model entries depend on 
free parameters for the particle potentials and level densi- 
ties, as well as for the 7 widths. These parameters have to 
be determined, tested, and improved by laboratory mea- 
surements. The HF predictions do also depend critically 
on the masses of the associated nuclei, whi ch are rather 
well known in the case of the p process l|Audi et al. I 
2003). A detailed description of HF techniques for de- 
termi ning nuclear reaction rate s are given in the litera- 
ture l|Rauscher fc Thielemannll2000lBoOll 1200^1 . The 
presently available HF models describe the astrophysi- 
cal reaction rates typically only within a factor of two. 
The specific limitations wit h respect to p-process anal - 
yses have been discussed bv lArnould fc Gorielv 1 12003). 
The observed discrepancies between experimental data 
and theoretical HF predictions seem to be largest in the 
case of a capture reactions which has been interpreted 
as the result of the use of insufficient a potential models. 
We therefore will give special consideration to the impact 
of uncertainties in (7, a) reactions on p-process reaction 
flow and p-nuclei abundance predictions. 

2. EXPERIMENTAL DATA 

Only a few of the required reaction rates have been 
investigated experimentally. While in most cases experi- 
mental and theoretical rates agree within the uncertain- 
ties d ue to the model parameters l)Arnould fc Oorielvl 
2003), significant discrepancies have been observed in 
some interesting cases. 

Neutron capture measurements fo r p-process studie s 
are available for part of the p nuclei iBao et al. 1 {2000), 
whereas methods for the determination of photodisintc- 
gration cross sections and reaction rates have been de- 
veloped only recent ly either using inverse Com pton scat- 
tered laser photo ns llUtsunom iya et al. 1200 ij) or electron 
bremsstrahlung i|Vogt et al. Il2001jh The latter method 
has been successfully applied to measure (7,n) photodis- 
integration react ions for the p pro cess including the p 
nuclei 190 - 192 Pt llVogt et al. ll20f)l . 19 M98 H& and thc 



proton magic nucleus 204 Pb ijSonnabend et al.l [20041. 
Good agreement with the HF predictions has been found 
in all these cases. Comparable data to test the reliability 
of (7,n) photodisintegration processes in the mass range 
below Z — 72 are not yet available, e.g. for nuclei near 
the N =50 and N =82 closed neutron shells. 

Considerable effort has been spent on determining the 
rates for (7,p) and (7, a) reactions via measurements of 
the inverse capture reactions. The approach has emerged 
as one of the major tools for testing the reliability of 
the HF predictions. The (p,7) rate for the p isotope 
96 Ru is on averag e a factor two below the HF results 
l)Bork et al. I Il998'l . On the other hand, experimental 
studies for the p nucleus 102 Pd indicate that the re- 
action rate for 102 Pd(p,7) 103 Ag is significantly higher 
than HF predictions l)Ozkan et al. I 120021k Recent p- 
nuclei related proton capture stu dies in the lower mas s 
range such as 74 - 76 Se(p /y'l 75 ' 77 B-r iGviirkv et al. 1^2 003) . 
or 92 < 94 Mo(p,7) 93 ' 95 Tc (ISaniter fc Kappeler 111 997^ 7^- 
pare rather well with HF predictions. Also recent proton 
capture reaction studies on neutron magic nuclei with 
N =50 such as 88 Sr fp.^ 89 Y iCalanopoulos et al. 112001 
and 89 Y(p, 7 ) 90 Zr ijTsagari et al. 112004^ seem to agree 
well with HF calculations within the parameter space 
of the model. In the range Z >50 experimental proton 
capture rates fo r p nuclei are not ava ilable except for 
112 Sn(p, 7 ) 113 Sb llChloupek et al. llT999h . which seems to 
be in excellent agreement with HF predictions, and for 
116 Sn(p,7) 117 Sb, which showed again significantly higher 
experimenta l values than predicted l|Ozkan et al. ll2002t 
IRapnll27)0l . 

Compared to neutron and proton associated capture 
and photodisintegration processes, larger deviations from 
HF predictions have been reported for a capture and 
photodisintegration reactions into the alpha channel. 
The proton and neutron binding energies for nuclei along 
the p-process path are typically very large, which war- 
rants the high level density conditions required for ap- 
plying thc HF model. For (7, a) reactions along the 
p-process path, however, the a binding energy is of- 
ten low and many nuclei are even a-unbound, spon- 
taneous a-decay being only suppressed by the respec- 
tive Coulomb barriers. In these cases, where devia- 
tions from HF predictions might be anticipated, thc 
number of experimental studies is small and limited 
to a capture and scattering reactions on N =82 and 
N =50 closed neutron shell nuclei along the p-process 
path. The measurement of the 144 Sm(a, 7) 158 Gd re- 
action (Q=-3.721 MeV) near th e N =82 closed neu- 
tron shell yielded a reaction rate ijSomoriai et al. 111998ft 
nearly o ne order of magnit ude lower than the HF pre- 
dictions l|Mohr et al. l lT997 h The calculations are based 
on a wide range of a potential parameters derived from 
extensive 144 Sm a elastic scattering dat a. Measure- 
ment s of a elastic scattering on 92 Mo l|Fiilop et al. I 
12001ft have been used to determine the 96 Ru(7, q) 92 Mo 
reacti on rate similar to the approach of iMohr et al. I 
l)1997ft . This prediction depends sensitively on the po- 
tential paramete rs and is considerably lower th an the 
HF calculations ijRauscher fc Thielemann! l200Cf) . Sim- 
ilar results have been obtained i n a capture mea- 
surements on the p nuclei 96 Ru l|Rapp et al. I 12002ft 
and 112 Sn l)Ozkan et al. 112002ft . where the experimen- 
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tal ra tes are also much lower than the HF predic- 
tions (Ra uscher fc Thielemann! 2000) for the tempera- 
ture range between 2.5 GK and 3.5 GK. Reasonably good 
agreement with theory was only found in the measure- 
ment of the (a, 7) cross section of 70 Ge ijFiilop et al. I 
1996). In additio n, recent (n,q) experiments on 147 Sm, 
143 Nd and 95 Mo llOledenov et al. 11200(1 IKoehler et all 
I200M IRapp et al. 1 12003ft have shown that present sta- 
tistical models overestimate the associated a-induced 
astrophysical rates by more than a factor of two. In 
the case of 147 Sm the a- widths distribution of the 
resonances showed in dications of non-statistical effects 
HKoehler et a,l. Il200l. 

In the following we want to study the impact of the 
nuclear reaction rates on the p-process yields within the 
framework of multi-mass zone simulations. After a short 
summary of the model parameters we will discuss the re- 
sults of the p-process nucleosynthesis simulations. These 
were performed by modifying the reaction rates within 
the uncertainty limits of the present experimental studies 
and of the theoretical predictions. The impact of these 
modifications will be analyzed in terms of the time inte- 
grated p-process reaction flux and in terms of abundance 
predictions for the p nuclei. A further goal of the present 
work is to identify the most critical reaction rates in the 
mass range A >57 with respect to the p abundances and 
to identify the reactions, which need more detailed ex- 
perimental investigation. 

3. MODEL PARAMETERS AND SIMULATIONS 

The present nuclear reaction network comprises more 
than 20000 reactions connecting about 1800 nuclei from 
hydrogen to bismuth. We have simulated the abun- 
dance evolution of the associated isotopes in the framc- 
work of a paramete rized type II SN shock front model 
ijRavet et al. 1119951) . The p process was investigated in 
14 different mass layers of the Ne/O burning zone of a 
25 Mq st ar for the temperature and density profiles pro- 
posed by l^shidaerar] (|2003) (see table [J). The se pro- 
files are comparable to the ones used bv ijRavet et al. I 
1995). It should be noted that a wide range of progen- 
itor masses are expected to contribute to the p-process. 
We chose here a 25 M Q star as a representative example, 
in part to facilitate compar i son with previous work, for 
examp l e |Gi>atiLiiLaL | (|2fJfJ0j) ; iWooslev. Heger fc Weaverl 
( 2002). R avet et al. I <|1995') surveyed the type II super- 
nova p-process for a wide range of models and find only 
a weak dependence of the final p-process on progenitor 
mass. 

Fig. |5]shows the temperature and density profile of the 
shock front passing through three of the layers reaching 
peak temperatures of 2.4 GK, 2.6 GK, and 2.96 GK, re- 
spectively. Table ^ lists the mass fraction of the star 
enclosed by each layer M„ as well as the peak tempera- 
ture and density reached in each of the layers. 

The final p-process abundances depend sensi- 
tively on the choice of the initial seed abundance 
l)Arnould & GorieTvl 12003ft . In the present approach 
we, therefore, kept the seed abundances fixed in order 
to study the sensitivity of the p abundances to nuclear 
reaction rates. The seed abundance distribution is 
determined by the s- and r-process history of the 
stellar material at the time of star formation, which 
is subsequently modified by in-situ nucleosynthesis 



during the evolution of the star considered. The latter 
aspect is particularly important since the Ne/O layers 
receive an abundance contribution in the mass region 
A «70-90 from the weak s process, which occurs 
during the preceding helium and carbon burning phases 
(|Kappeler et al. lll994l:lThe et al. 112000ft . This s-process 
component depends on the mass of the star but also on 
the neutron irradiation provided by the 22 Ne(a,n) neu- 
tron source. There has been some speculation that the 
22 Ne(a,n) so urce might be more e fficient than previously 
anticipated l)Costa et al. I 12000ft . which would cause 
a substantial enhancement in the abundances of the 
p-only nuclei 92 Mo and 96 Ru . A significantly enhanced 
22 Ne(a,n) reaction rate would, however, result in an 
enormous overproduction of s-process nuclei, in sharp 
confli ct with the ob served galactic abundance patterns 
ijHeger et al. 1 1200*211 . Also, new measur ements of the 
22 Ne (a,n) cross section at low energies ijJaeeer et al. I 
12001ft and a de tailed analysis of the 25 Mg(n,7) 26 Mg 
reaction channel ijKoehler 112002ft provide stringent limits 
for th e uncertainty of this reaction rate l)Karakas et al. I 
12004ft and exclude the suggested enhancement factor. 
The seed abundances for the present study (Fig. |3J) 
were adopted from a pre-supernova e v olution model 
of a 25M^ star jRavet et al. I 119951 iRavet I l2002t 
lArnould fc Gorielv 1 12003ft to facilitate comparison with 
previous work. 

For modeling the p-process nucleosynthesis yields 
all n-, p-, and a-induced capture rates as well 
as their inverse photodisintegration rates were based 
on the HF predictions of the NON-SMOKER code 
ijRauscher fc Thielemann! 12000ft . In addition, re- 
action rates on light n uclei with Z <8 were in - 
cluded from the work o f iCaughlan fc Fowler I j!988ft . 
[ Wiescher et al. I lll989ft . iRauscher et al. I ()1994ft . and 
iHerndl et al. I l)1999ft to account for the impact of light 
particle capture reactions on the proton, neutron, and al- 
pha budget. The abundances of the following 35 p nuclei 
were investigated in the calculations: 74 Se, 78 Kr, 84 Sr, 
92 Mo, 94 Mo, 96 Ru, 98 Ru, 102 Pd, 106 Cd, 108 Cd, 112 Sn, 
113 In, 114 Sn, 115 Sn, 120 Te, 124 Xe, 126 Xe, 130 Ba, 132 Ba, 
136 Ce, 138 La, 138 Ce, 144 Sm, 152 Gd, 156 Dy, 158 Dy, 162 Er, 
164 Er, 168 Yb, 174 Hf, 180 Ta, 180 W, 184 0s 190 Pt, and 
196 Hg. Following the method described in IRavet et al. I 
( 1995) the mass fraction X^„ of each p nucleus i in each 
p-process layer n was calculated independently for all 14 
layers. The total mass m-i of nucleus i was determined by 
the respective contributions weighted by the mass of each 
zone delimited by its two neighboring p-process layers, 

n— 13 

n>l 

with M„ representing the mass of the star within layer 
n. 

3.1. p-process abundances 

The efficiency of a particular nucleosynthesis process 
in contributing to the observed solar abundance distribu- 
tion can be expressed by the overproduction factor (i^), 
which compares the produced abundance of a given iso- 
tope i with the observed solar abundance. For the anal- 
ysis of the present calculations an overproduction factor 
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has been denned as 



(Fi) = 



(2) 



with the total mass M tot — J2 n >i M n — M n -i = Mi 3 — 
Mo obtained by the sum over all p-process layers . The 
solar abundance ma ss fractions were taken from 
lAnders et al~lft989|) . The averaged overproduction fac- 
tor for the 35 p-only nuclei is 

(Fi) 



Fn = 



^ 3.1 



(3) 



The normalized overproduction factor (Fi)/F is by def- 
inition equal to unity when the simulate d abundances 
matc h with the observed solar values l)Ravet et al. I 
1995). 

Fig. 0] shows the normalized overproduction (Fi) /Fq 
as a function of mass number A. Given the uncertainties 
of the present model and the uncertainties of the nuclear 
input the calculated p-abundance pattern compares rea- 
sonably well with the observed solar abundances. The 
comparison of the results displayed in Fig. 0] with previ- 
ous work relies on an identical seed abundance distribu- 
tion and on similar statistical model based reaction rates. 
Within the expected model specific discrepancies our cal- 
culations yield a p-nuclei abundance pattern, which is 
very similar to the previous res ults. Our prediction for 
190 Pt disagrees with the result of Rav eTet al. Illl995ft but 
is con sistent with the value given bv lArnould fe GorielvH 
( 2003). This may reflect the use of different reaction rate 
compilations. 

The observed abundances of the p nuclei could be re- 
produced in most cases within a factor of three. Par- 
ticular exceptions are the notorious underproduction of 
92 - 94 Mo and 96 ' 98 Ru and the deficiencies of U3 In, 115 Sn, 
and 138 La. Also the much debated 180 Ta abundance is 
found in agreement with previous calculations. How- 
ever, this case has been treated without distinguishing 
between ground and isomeric state. Accordingly, the 
present value represents an upper limit for the p contribu- 
tion to that isotope, still compatible wit h the significant 
s-prq cess contribution from AGB stars ijKanueler et al.1 
2004). 

The low production of the light p nuclei 92,94 Mo and 
96,98 Ru has been observed in all previous p-process stud- 
ies lIRavet et ah 1 179951 lAmonld fc Gorielv I Eool and 
remains as the main enigma in p-process simulations. 
Possible additional nucleosynthesis scenarios for produc- 
ing these isotopes have been discussed and include a 
i/-induced abundance componen t originating from th e 
neutrino-driven supernova shock ijHoffman et al. l fl996). 
the production in the neutron-rich, alpha-ric h freeze 
out n ear the mass cut in type II supernovae l)Mever I 
12003ft . and finally the independent production mecha- 
nism by the rp p rocess in the outer l ayers of accret- 
ing n eutron stars l)Schatz et al. I Il998t iDaunhas et al. I 
2003). Within the framework of the p-process models 
possible explanations have been relat ed to modifications 
in the seed abundance di stribution ((Costa et al. 1 1200(1 
lArnould fc Goriefvl 12003ft . as discussed before, or were 
suspected to result from the nuclear physics parameters 
used for calculating the associated reaction rates. 

Also the underproduction of 113 In, 115 Sn, 138 La, 
152 Gd, and 164 Er has been observed in previous calcu - 
lations l|Ravet et al. 1119951 lArnould fc Gorielv 1 12003|) . 



No other nucleosynthesis source has been identified 
for 113 In since it i s effectively shielded from the s- 
and r-process path ijNemeth et al. 1 119941 iTheis et al. I 
11998(1 . For the underproduced 115 Sn the s process 
provides an additional 5% production through the se- 
quence 114 Cd(n,7) 115 Cd m (/r) 115 In m feeding the 1/2" 
isomeric state in 115 In™ 1 , which su bsequently decays 
through a weak 0~ branch to 115 Sn l|Beer et al. 1119891 
iNemeth et al. lll994ft . More significant contributions to 
113 In and 115 Sn (12% and 43%, respectively) are ex- 
pected from the r process via the /3-unstable isomers 
in 113 Cd and 115 In, which are populated in the post-r- 
proce ss decay chains ((Nemeth et al. lfl994t ITheis et al. I 
1998). Nevertheless, the origin of the rare In and Sn 
isotopes represents an unsolved puzzle. 

The abundances of 152 Gd and of 164 Er, however, 
are understood to result mostly from s-process nucle- 
osynthesis. In total 71% of 152 Gd are produced by 
the s -process branching at 151 Sm ijAbbondanno et al. I 
12004ft . The here observed abundances for 152 Gd orig- 
inate mainly at relatively cool temperatures Tg < 2.0 
typical for the outer p-process layers. Similarly, 164 Er 
is underproduced by the p process but more than 90% 
can be explained by feeding through the high temper- 
ature s-process branching via bound-state /3-decay of 
163 Dy to 163 Ho with the subsequen t neutron capture 
sequence 163 Ho(n ^) 164 Ho(/r ) 164 Er llaag fc Kanneler I 
199$ IBest et al. ll200l . 

The nucleus 138 La is systematically underpro- 
duced in all supernov a based p-process models 
l|Arnould fc Gorielv 1 12003ft . Alternative models rela ted 
to energetic stellar pa rticles llHaineb ach et a l. Ill976j) or 
neutrino interactions <|Wooslev et al. 1 11990ft have been 
developed to explain the observed 138 La abundances. A 
de tailed analysis o f the 1 38 La problem has been reported 
bv iGorielv et al.1 (|200 If) . In terms of nuclear reaction 
rates it was shown that plausible estimates for the re- 
lated uncertainties can not explain the underproduction 
of 138 La. This uncertainty could be removed, however, 
if the theoretical production rate 139 La(7,n) and the de- 
struction rate 138 La(7,n) could be replaced by experi- 
mental data. 

The final p abundances depend on the p-process reac- 
tion flux and reaction branchings, which determine the 
feeding and depletion of the various p nuclides. In the fol- 
lowing section we discuss the characteristic flux patterns 
as a function of mass layer and/or peak temperature. 

3.2. p-process flux 

The time integrated reaction flux per mass layer pro- 
vides information about the main reaction path during 
the nucleosynthesis event and serves as a tool for moni- 
toring the effects of nuclear structure parameters such as 
shell closure or deformation on reaction path and reac- 
tion branchings. The time integrated reaction flux cor- 
responds to the net number of reactions between two 
isotopes i and j integrated over a certain period of time 
At = <i — to. This flux is defined by 



fid 



it, 



dt (i^j) 



dt (j^i) 



dt 



(4) 



with the isotopic abundances Yi — Xi/Ai (mass fraction 
divided by mass number). In the following we refer to 
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the time integrated reaction flux simply as reaction flux. 
The maximum reaction flux fij defines the main reaction 
path along which nucleosynthesis will take place. 

The reaction flux has been calculated separately for 
each investigated layer and depends strongly on the re- 
spective peak temperatures. The temperature and den- 
sity profiles are shown in Fig. [21 Fig- El shows the reac- 
tion flux integrated over a 1 s time interval covering the 
entire shock-driven temperature peak of T g =2.96 shown 
in Fig. [21 During this period the seed nuclei in the range 
62< Z <83 are processed by (7,n) reactions towards the 
neutron deficient side of stability (upper part). With 
increasing neutron binding energy of the reaction prod- 
ucts (Fig. 0) the abundance distribution is driven to- 
wards lower masses by the competing (7,01) channel. Be- 
cause the neutron deficient isotopes above N =82 are 
a-unbound as shown in Fig. [7| spontaneous decay is 
only inhibited by the high Coulomb barriers. For neu- 
tron magic nuclei the rapid increase in neutron and a 
binding energy reduces the (7,11) and (7, a) photodisin- 
tegration rates significantly and forces the reaction flux 
into the (7,p) channel along the N =82 closed shell nuclei 
towards lower masses. In particular the p-process iso- 
tope 144 Sm is formed by feeding through the 147 Eu(7,p) 
and the 148 Gd(7, a) channels, but is depleted through 
144 Sm(7,n) photodisintegration. 

The middle part of Fig.EJshows the integrated reaction 
flux for 48< Z <62. While in this mass range the overall 
reaction path is still driven by (7,11) reactions towards 
the neutron deficient side of stability, the flux pattern 
is characterized by a strong (7,p) reaction component 
towards stability for N <82 since the rapid increase in a 
binding energy (Fig. EJ inhibits the (7,0;) channel. The 
(7,0;) processes regain a more competitive role for nuclei 
with Z <58 and provide a strong reaction flux towards 
the even-even proton magic Sn isotopes with Z =50. The 
main reaction channels pass through 112 Sn and 114 Sn, 
which are both fed via (7,11) and (7, a) reactions. 

The p-process flux below Z =50 (lower part of Fig. [5j 
is still characterized by (7, a) reactions feeding the p nu- 
clei 106 ' 108 Cd, 102 Pd and eventually 96 Ru and 92 Mo. At 
N =50 the (7, a) flux is largely diminished because the 
increasing a binding energy efficiently reduces the (7, a) 
channel. An additional strong (7,11) and (7,p) flux oc- 
curs towards N =50 closed neutron shell nuclei, which is 
followed by (7,p) photodisintegration along the N =50 
isotone chain. Compared to the integrated flux in the 
region Z >50 the overall flux is significantly reduced be- 
cause of the high single particle binding energies in this 
region near the Z, N =50 closed shells. The figure indi- 
cates that the abundances of the underproduced p nu- 
clei such as 92 Mo and 96 Ru are mainly determined by 
the (7,n) reaction flux feeding and the (7,p) flux deplet- 
ing these isotopes. In addition there is a significant flux 
through the 96 Ru(7, a) 92 Mo reaction. The final abun- 
dances of the p nuclei 94 Mo and 98 Ru depend mainly on 
the (7,n) rates depleting these isotopes with 94 Mo being 
mainly produced at the lower temperature mass zones. 
Only a few (7, a) rates contribute to the reaction flux 
at Z <50, mainly due to the rapid increase in a-binding 
energy as shown in Fig. 

In the mass range below the N =50 neutron shell the 
nuclei become more and more resistant against photo- 
disintegration because of their high neutron and proton 



binding energies (Figs. E]and|S|. The reaction flux is 
characterized by a complex pattern of single neutron, and 
proton capture and their inverse dissociation reactions. 
These processes are accompanied by (p,n) and inverse 
(n,p) reactions because of their typically low thresholds 
of w 0.5 - 1.0 MeV. The reaction path remains close to 
stability. 

A case of p-process nucleosynthesis in a layer reaching 
only a peak temperature of Tg=2.44 is displayed in Fig. 
El For isotopes with 62< Z <82, shown in the upper part 
of the figure, the reaction flux is dominated by (7,11) reac- 
tions, which drive the initial seed abundance distribution 
towards the neutron deficient side of the line of stability. 
After the shock front passed, the produced radioactive 
isotopes decay back to the line of stability. Only minor 
changes are anticipated for the overall abundance distri- 
bution in this mass range. Besides the (7,n) reactions a 
few (7, a) reactions can be observed processing neutron 
deficient isotopes in the Hf to Pt range towards lower 
masses. At these temperatures the production of 162 Er 
and 184 0s (Fig. llOfl is limited to the (7,11) reactions 
along the Z =68 and Z =76 isotope chains, respectively. 
In the case of 164 Er, no feeding through 168 Yb(7, a) 164 Er 
seems to be taking place. The p contribution to 152 Gd 
depends on the (7,n) reactions along the Z =64 isotope 
chain, but is further depleted through 152 Gd(7, a) 148 Sm. 
The N =82 closed shell isotope 144 Sm is mainly fed by 
148 Gd(7, a) 144 Sm reactions. 

In this mass layer the time integrated flux between the 
N =82 and Z =50 closed shells is again dominated by 
(7,n) reactions subsequently balanced by inverse (n,7) 
neutron capture reactions as shown in Fig. E] There- 
fore, the abundance distribution changes only along the 
isotopic chain. Both (7,p) and (7, a) reactions are negli- 
gible. The p nuclei in this mass range are not affected by 
the reaction flux with two exceptions, 138 La and 138 Ce. 
This pattern continues also towards lower masses down 
to the N =50 neutron closed shell. The (7,n) reaction 
flux associated with the 92,94 Mo and 96 ' 98 Ru nuclei is 
negligible because of the high neutron binding energies. 
These results suggest that the p-process production of 
these isotopes is confined to higher temperature layers. 

For nuclei below the N =50 closed neutron shell the re- 
action pattern differs distinctively from the higher mass 
regions. Similar to the flux at higher temperatures, pho- 
todisintegration into the neutron and proton channel and 
their inverse capture processes as well as (p,n) and in- 
verse (n,p) reactions define the reaction path and deter- 
mine the final abundance distribution in the mass range 
below N =50. 

As shown in Figs. El and El the integrated p-process 
reaction flux is considerably different between high 
(Tg=2.96) and low temperature (Tg=2.4) layers. In par- 
ticular the nucleosynthesis of higher mass p nuclei such as 
156 Dy, 162 Er, 164 Er, 168 Yb, 174 Hf, 180 W, 184 0s, and 190 Pt 
is temperature sensitive because of the strong tempera- 
ture dependence in the (7, a) reaction rates, which deter- 
mine the reaction flux pattern. The critical peak temper- 
ature for these reactions is Tg=2.6, where the reaction 
path is still largely characterized by a strong (7, a) flux 
towards lower mass isotopes. Below Tg=2.6 the (7, a) 
reaction rates are typically too weak to warrant rapid 
depletion of high Z isotopes. This causes an enrichment 
of the p nuclei through (7,11) feeding from the initial seed 
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abundance distribution. For temperatures above 2.6 GK 
the abundances of the heavy p nuclei are reduced due to 
significant (7, a) depletion, which increases exponentially 
with temperature. 

The reaction flux depends on the initial seed abun- 
dance as well as on the reac tion rates. It has been argued 
ijArnould fc Gorielv"ll2003|) that the reaction rates have 
only a limited influence on the final p-process abundance 
distribution and in particular will not solve the mystery 
of the light p-nuclei abundances. In the following we 
seek to investigate in an independent systematic study 
the direct impact of p-process reaction rates. 

4. THE INFLUENCE OF REACTION RATES 

The influence of reaction rates on the final p-process 
abundance distribution was studied first by introduc- 
ing global enhancement and reduction factors for dif- 
ferent reaction channels and by probing the sensitivity 
of the different (7,n), (7, a), and (7,p) photodisintegra- 
tion branches. All other input parameters and condi- 
tions remain unchanged in the model. The approach of 
collectively changing all rates helps to identify the mass 
range where rates have a direct impact on the p-process 
flux and p-process abundance predictions. The final p- 
abundances are sensitive to reaction rates because they 
determine the balance of feeding and depletion flows for 
each p-nucleus. In addition, changes in reaction rates 
change branching ratios when several reaction channels 
compete, for example (7,11) with (7,p) in the element 
range up to Z = 64 (see Fig. 5). In addition, such 
branchings are sensitive to nuclear structure effects such 
as shell closures. 

4.1. Collective change of neutron- proton- and a-rates 

As a first test of the sensitivity of the p-abundance dis- 
tribution the reaction rates for all neutron-induced pro- 
cesses and their inverse photodisintegration reactions on 
nuclei with A >57 were enhanced or reduced by a factor 
of three. In each case the p abundances were calculated 
and the overproduction factors were compared with the 
results based on the standard set of HF rates used in this 
study l|Rauscher &: Thielemann ll2000j) . Sensitivity stud- 
ies for proton- and a-induced rates were performed in a 
similar way. Fig. EUshows the respective overproduction 
factors calculated with the modified rates in comparison 
with the ones obtained with the standard HF rates. 

The impact of the collective change of all (n,7) and 
of the inverse (7,11) rates on the ^-nuclei abundances is 
shown in the top part of Fig. ^| One observes a corre- 
lation over the entire mass range between the rates for 
the neutron channel and the p-process yields. This is 
not too surprising because (7,n)-reactions are active in 
all of the investigated layers where the p process is taking 
place. The overall abundances of the p nuclei change by 
less than a factor of two on average. Since the reaction 
flux below N =50 is driven by capture and photodisinte- 
gration reactions the abundances of p nuclei with A <80 
correlate more strongly with changes of the reaction rate. 
A more pronounced correlation can also be observed for 
162 Er, 164 Er, and 168 Yb, which are sensitive to the (7,n) 
photodisintegration rates in the mass layers exposed to 
lower peak temperatures 2.76 < Tg < 2.32 as indicated 
for the example of 162 Er in Fig. ^2 

The middle part of Fig. ^] displays the sensitivity of 



the predictions to global changes in the proton-induced 
and inverse reaction rates. No impact can be observed in 
the mass range A >110 despite the fact that the reaction 
flux pattern shows a pronounced (7,p) flux component 
between N =50 and N =82. In our model the only ex- 
ception is 126 Xe, which is produced by /3-decay of 126 Ba. 
The final abundance depends strongly on the (7,p) pro- 
duction chain of 126 Ba in the higher temperature mass 
layers. A more pronounced effect is obtained in the lower 
mass range A <110 where the p nuclei with Z <50 and 
N >50, such as 92 Mo and 96 Ru, are mostly depleted by 
(7,p) photodisintegration (Figs. ElandEJ). A reduction 
in rate is therefore directly correlated with an increase in 
abundance. The high neutron and alpha binding ener- 
gies of the neutron magic nuclei with N = 50 terminate 
the (7, n)-reaction flux at 92 Mo (see Figs. ElandEj). At 
this point the flux continues only via the (7,2?) rate be- 
cause of the relatively low proton binding energy (see 
Fig. EJ. Hence, the 92 Mo abundance is strongly deter- 
mined by the rate of the 92 Mo(j,p) reaction. In contrast, 
the 94 Mo(7,p) reaction has no effect on the 94 Mo abun- 
dance since the much lower neutron binding energy of 
this isotope results in a dominance of the (7, n) chan- 
nel. Similarly, 96 Ru and 98 Ru are both strongly bound 
against (7, n) reactions (see Fig. |SJ, but 96 Ru has a suf- 
ficiently small proton binding energy as well to exhibit 
a clear sensitivity to changes of the ("f,p) rate. There- 
fore, the final abundances of the p nuclei 94 Mo and 98 Ru 
remain unaffected by changes in the (7,p) rates. The 
abundances of p nuclei in the low mass range, such as 
74 Se and 78 Kr correlate inversely with the strength of 
the proton capture reactions (chapter 4.4). 

Global variation of the (a, 7) and (7,0;) rates has a 
strong impact on the p-nuclei abundances above A =140 
and N >82 as shown in the bottom panel of Fig. 1101 
In this range the reaction pattern is strongly affected 
by (7, a) photodisintegration reactions (Figs. 03and|nj|. 
With increased rates the higher mass p nuclei are by- 
passed and the material is processed towards the lower 
mass p nuclei such as 144 Sm, 156 Dy and 162 Er. If the re- 
action rates are lower, processing towards these p nuclei 
is less efficient and causes a relative enrichment in the 
higher mass nuclei such as 174 Hf, 180 W, and 190 Pt. The 
abundance of 152 Gd remains unaffected by changes in the 
a capture or emission rates. In the case of 164 Er, varia- 
tion of the (7, a) rates leads always to a decrease in abun- 
dance. If the rate is reduced the direct feeding through 
168 Yb(7, a) 164 Er is suppressed, and if it is increased the 
abundances of higher mass feeding isotopes are reduced 
by faster processing towards lower mass N =82 isotones. 
In the lower mass range the abundance of 96 Ru is by far 
the most sensitive to changes in the (7, a) rates. Since 
the depletion depends significantly on the 96 Ru(7, a) 
rate, a corresponding enhancement of this rate causes 
a substantial reduction in the 96 Ru abundance. This is 
similar to the case of enhancing the 96 Ru(7,p) depletion 
rate discussed above. 

As pointed out before, experimental studies indicate 
deviations of up to one order of magnitude with respect 
to HF predictions for a capture and their inverse photo- 
disintegration processes. We therefore performed a sec- 
ond p-process simulation with the (7, a) rates modified 
by factors of 0.1 and 10, respectively, to account for a 
broader range of uncertainty. Fig. El shows that these 
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changes result in a similar pattern as observed in the 
previous study (bottom panel of Fig. 11(1)1 . A signifi- 
cant change can only be discerned for 156 Dy, 162 Er, and 
to a lesser extent for 190 Pt. The increase with higher 
reaction rates is due to enhanced feeding of 156 Dy and 
162 Er through (7, a) reactions in the mass zones with 
lower peak temperatures (Fig. ITT|) . 

We have focused so far on global changes of reaction 
rates. In the following we discuss the impact of reac- 
tions feeding or depleting only the p nuclei in order to 
determine how critical these reactions are for the overall 
abundance predictions. 

4.2. Impact of rates for p-nuclei reactions 

Similar to the discussion before, the rates directly feed- 
ing or depleting the p nuclei were changed by factors of 
three to see if the calculated p-process abundances are 
determined by the global reaction flux or if they are as- 
sociated with these specific reactions. This will have con- 
sequences for the identification of reactions to be selected 
for further experimental studies. 

Changing the neutron capture rates and the respec- 
tive photodissociation rates into the neutron channel by 
a factor of three produced only negligible changes in the 
resulting p abundances compared with the results of the 
global rate changes. In particular the change of (n,7)- 
rates on the p nuclei showed almost no influence while 
the change of the (7,n)-reactions showed mainly a slight 
sensitivity between 80< A <160. Since (7,n) photodisin- 
tegration processes provide the main feeding path for the 
transformation of the initial seed nuclei towards neutron 
deficient isotopes no specific outstanding reactions could 
be identified in this analysis. The reason is that the 
p-process reaction flow proceeds through neutron defi- 
cient nuclei several mass units away from stability. With 
increasing neutron separation energy the (7,n) reaction 
rates decline rapidly and the reaction flow is carried by 
(7,p) and (7, a) channels. It is important to identify these 
branching points which determine the overall p-process 
reaction flow pattern towards lower masses. The branch- 
ing points are clearly temperature dependent and differ 
between the different burning zones. The highest sensi- 
tivity appears to be correlated with (7,n) reactions near 
the N = 50, 82 closed neutron shells where the (7,n) 
reaction flow changes into a (7,p) dominated reaction 
flow pattern. In addition, the weak sensitivity of the p 
abundances to the (7,11) reactions feeding the p nuclei is 
not unexpected. These reactions typically occur on odd 
N nuclei with lower Q-values and have therefore higher 
rates compared to neighboring (7,11) reactions on even 
N nuclei. They are therefore not expected to be partic- 
ularly important bottle-necks in the p-process reaction 
flow. 

In the case of (7,p) dissociation reactions on p nuclei 
and their inverse processes, a modification of these rates 
has a direct impact on the abundances of p nuclei in the 
mass range A <100. Fig. 1131 compares the overproduc- 
tion factors of p nuclei based on the modified rates to the 
ones based on standard HF rates by showing the ratio as 
a function of atomic mass number A. The squares (factor 
1/3) and crosses (factor 3) reflect the global changes of 
all proton related reaction rates discussed before (see also 
the mid part of Fig. I1U|) . whereas the solid and dashed 
lines denote the corresponding ratios obtained by chang- 



ing only the (7,p) dissociation rates on p nuclei and their 
inverse (p,7) reaction rates. The results indicate that the 
reactions directly associated with p nuclei are the most 
critical ones, though in a few cases (A = 74, 102, 106, 
108) reactions on non-p nuclei are also important (see 
also table 2). The figure clearly shows that the (7,p)- 
rates determine significantly the p abundances of 74 Se, 
78 Kr, 84 Sr, 92 Mo, and 96 Ru since the overabundance ra- 
tios scale inversely with the rate scaling factors. Proton 
capture rates have almost no influence on the p abun- 
dances, they play only a role for the light p nuclei 74 Se 
and 84 Sr. The comparison of the overabundance predic- 
tions based on the globally changed rates and selectively 
changed rates suggests that for reactions involving pro- 
tons the individual feeding and depleting processes of p 
nuclei are significant. Fig. 1131 indicates that the abun- 
dances of 78 Kr, 92 Mo, and 96 Ru depend sensitively on the 
photodisintegration of these nuclei. It has been pointed 
out before that particularly the 92 Mo(7,p) 91 Nb rate de- 
termines the 92 Mo / 92 Nb abundance rat io predicted by 
supernova models (Da uphas et al. 1 12003). 

Also in the case of a capture and a disintegration re- 
actions the results indicate that the produced overabun- 
dances are closely correlated to the individual feeding 
and depleting reactions of the light p nuclei. Correla- 
tions in specific cases have already been found in Fig. 1121 
with ( 7 ,a) reactions depleting 74 Se, 96 Ru, 120 Te, 122 Xe 
and ( 7 ,a) reactions feeding 102 Pd, 106 Cd, 144 Sm, 156 Dy 
and 162 Er as discussed in the chapter 4.1. 

The present simulations yield low abundances for the 
Mo and Ru p isotopes, in agreement with previous cal- 
culations. The difficulty in solving this problem within 
the astrophysical model and seed distribution discussed 
here comes from the fact that 92 ' 94 Mo and 96 ' 98 Ru make 
up about 43% of the total solar p nucleus abundances. 
Therefore, increasing their overproduction significantly 
requires a very efficient conversion of a significant frac- 
tion of all the heavier seed nuclei into Mo and Ru. As 
Fig. ITT1 shows, an absolute overproduction factor of the 
order of 10 — 90 would be needed for the Ru and Mo p 
isotopes to be in line with heavier p nuclei, assuming that 
specific isotopes are roughly produced in layers of similar 
mass. For comparison, conversion of the entire Z > 42 
seed into 92 ' 94 Mo and 96 ' 98 Ru would be required to bring 
their overproduction factors up to an average of 70. Even 
if 92 ' 94 Mo would be produced in separate layers without 
coproduction of 96,98 Ru the maximum achievable over- 
production factor would onl y be around 110. I n con- 
trast to the earlier findings of Ravet et al. I l)1990j) using 
a parametrized p— process model this would be in prin- 
ciple sufficient, but it is difficult to see how such a major 
increase in Mo and Ru production efficiency could come 
about. Nevertheless, to investigate this issue further, we 
decreased the 92 - 94 Mo(7,p) and 96 > 98 Ru(7,p) rates to 10% 
of their HF prediction. These reaction rates are the main 
destruction mechanism for 92,94 Mo and 96,98 Ru. The re- 
sulting change in abundance is shown in Fig. 1141 Only a 
moderate enhancement of the associated abundances is 
observed. Accordingly, reaction rates are not responsible 
for the underproduction of the light p nuclei. 

4.3. Critical p-process reaction rates 

The simulations discussed in the previous sections indi- 
cate that within the present model relatively few selected 
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reaction rates have a major impact on the final p abun- 
dances. The critical rates are typically associated with 
a strong feeding flux or with a particular branching in 
the reaction flux feeding or bypassing the p nuclei (Fig. 
EJ). The (7,n) reactions control the overall feeding of the 
p-process flux from the weak s-process seed distribution, 
thus affecting the final j»-process abundance distribution 
over the entire mass range. This influence is particularly 
visible in the range above N =82 as suggested by Fig. ITCH 
(top panel). In the following we want to concentrate on 
the discussion of the (7,p), (n,p), and (7, a) branchings 
which divert the p-process flux towards lower masses. 
We identified a set of (7,p) and (n,p) reactions and the 
respective inverse processes to investigate their specific 
impact on the simulations. Table [5] also includes those 
reactions in the mass range 70< A <78, which seem to 
affect the abundances of the very light p nuclei below 
92 Mo. 

A change in the rates by a constant factor can result 
in opposite sensitivities as can be seen for 74 Se and 78 Kr 
in the middle panel of Fig. 1101 For proton induced re- 
actions and their inverse reactions the 78 Kr abundance 
depends only on the depletion reaction 78 Kr(7,p) (see 
top part of Fig. I14|) . hence an increase of the proton- 
induced rates results always in a lower 78 Kr abundance. 
In contrast, the 74 Se abundance depends not only on 
the 74 Se(7,p) rate since it is produced in the high tem- 
perature layers where the flow via 74 Ge(p, 7) 75 As and 
75 As(p, n) 75 Se(7, n) leads to an enhancement of the 74 Se 
abundance as these rates are increased. 

The influence of the selected (7,p) and (p,n) reactions 
as well as their inverse processes is demonstrated in the 
top panel of Fig. 1151 where the abundance modifications 
resulting from the global change of all (7,p) and (n,p) 
reaction rates and their inverse processes by a factor of 
three (mid panel of Fig. I10J) are compared with the ones 
based on the change of the rates listed in Table [3 The 
solid and dotted lines indicate the consequences of in- 
creasing and decreasing these rates by a factor of three, 
respectively. Within the general uncertainties the ob- 
served abundance pattern agrees very well with the abun- 
dance pattern resulting from a global change of all rates 
within the same boundaries and shows that the impact of 
(7>p) photodisintegration processes is confined to the p- 
process flux in the range Z <50 with the only exception 
of 126 Ba(7,p) 125 Cs. The impact of the (7,p) reactions on 
the abundance predictions for the light p nuclei 92 Mo and 
96 Ru is substantial but - within the chosen boundaries - 
not strong enough to explain the large under-production 
of these crucial p nuclei. Sensitivities to the uncertain- 
ties in the reaction rates of the (p,n) processes are mainly 
observed in the p-process range N <50 where photodis- 
integration plays a less dominant role than in the higher 
mass regions. 

The high seed abundances in the mass region from 
Hg to Pb (Fig. |3J) are efficiently shifted by (7, a) reac- 
tions towards lower masses, particularly at temperatures 
of Tg = 2.6. This mass flow is considerably enhanced 
by increasing the rates of the relevant reactions listed 
in table 3. For example, the abundances of the p nu- 
clei 156 Dy and 162 Er exhibit a strong sensitivity on the 
160 Er( 7 , a) 156 Dy and 166 YbJ/7, a) 162 Er reactions, respec- 
tively. As shown in Fig. 1141 increasing and decreasing 
these rates by factors of three results in a correspond- 



ing decrease and increase of the overproduction factors 
of 156 Dy and 162 Er. However, the global change of the 
(7, a) rates by a factor of ten as illustrated in Fig. ^] 
leads to an inverse trend in the overproduction factors 
for 156 Dy and 162 Er, which are strongly enhanced also 
if the rates are increased. This behavior reflects the en- 
hanced mass flow from the abundant seeds in the Hg to 
Pb region. In summary, the importance of the (7, a) 
reactions calls for more systematic studies in the mass 
region A > 140. 

Table lists the (7, a) reactions, which carry a sub- 
stantial reaction flux in the p process (see Fig. EJ . Again, 
modifying these reactions within the given uncertainty 
range causes significant changes in the abundances of the 
p nuclei. This is demonstrated in the lower panel of Fig. 
1151 bv comparing the abundance modifications resulting 
from the global change of all (7, a) reactions in the net- 
work by a factor of three (bottom panel of Fig. I10[l with 
the ones based only on the modified rates listed in Table 
01 Again, the solid and dotted lines indicate the response 
to an increase and the reduction of these rates by a factor 
of three, respectively. 

The changes in p-nuclei abundances based on the mod- 
ification of the selected rates in Table El agree well with 
the global change of all (7, a) rates. The figure sup- 
ports the argument that the impact of the (7, a) pho- 
todisintegration processes is mainly visible in the higher 
mass p-process range above N =82. This underlines the 
particular importance of (7, a) rates for the mass range 
N >82 except for the case of the light p nuclei 96 Ru and 
7 Se, where the predicted abundance is correlated with 
the depletion via 96 Ru(7, a) 92 Mo and 74 Se(7,a) 70 Ge re- 
spectively, as indicated in Fig. El 

The analysis of the abundances for p nuclei with N >82 
indicates that the low abundances of the p nuclei 156 Dy, 
162 Er, and 190 Pt might be due to the (7, a) rates in that 
mass range, while p nuclei such as 152 Gd, 158 Dy, and 
1 Er depend more on the strength of the particular (7,11) 
feeding from the seed abundance distribution. 

The results shown in Fig. ED suggest that compa- 
rably few rates involving charged particles have a di- 
rect impact on the abundances of the p nuclei. These 
critical rates, which are summarized in Tables and 
El carry most of the charged particle related uncertain- 
ties for simulations within the discussed p-process model 
frame. Since all of the se rates are presently based on 
global HF calculations l)R auscher fc Thielemann I [20001 
I2001L 120041) . experimental confirmation is necessary for 
reducing the inherent nuclear model uncertainties. This 
serves not only the purpose of obtaining better data, but 
since p-process abundances are sensitive to changes in 
these specific rates, a better knowledge of those rates 
will help to improve the fine-tuning for modeling of the 
supernova shock front traversing the O-Ne shell of the 
pre-supernova star. 

Our list of potentially i mportant reaction s largely dif- 
fers from the one given bv lRauscherl l|2006[K This is not 
surprising because the latter work identifies reaction flow 
branchings in a more simplified approach and does not 
follow the hydrodynamic evolution of different p-process 
layers. In addition, we have not limited our analysis to 
branchings, and we have also included other types of re- 
actions such as (n,p) or (p,n). 

There are experimental data for some of the listed (p,n) 
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reactions. Fig. 1161 shows the experimental data for the 
75 As(p,n) and the 85 Rb(p,n) reaction in comparison with 
Hauser Feshbach predictions. Only three data points 
have been determi ned for 75 As(p,n) in the energy range 
from 3 to 5 MeV ijKailas et al. Ill979ft . The experimen- 
tal results are in excellent agre ement with the statistical 
Haus er Feshbach predictions ijRauscher fc Thielemann I 
1200 1|) . On the other hand, the 75 As(p,n) cross sections 
obtained by a measurement of 75 As(p,xn) reactions in 
the energy range of 3 - 45 MeV lIMushtaaet et al. lll98l 
show significant discrepancies to the model predictions 
for the astrophysical low energy range. Experimental 
cross section data for the 85 Rb(p,n) reaction are avail- 
able for a wide energy ran ge from 3 MeV up to 100 
MeV ijKastleier et al. 1 12002'). The experimental results 
are o n average lower than the Hause r Feshbach predic- 
tions ijRauscher fc Thielemann ll2001|k on the other side, 
the data are handicapped by huge experimental uncer- 
tainties. In both cases independent experimental verifi- 
cation of the cross sections for the astrophysical energy 
range is desirable. 

Different experimental techniques will be necessary for 
a complete study of these critical reactions. As men- 
tioned before, in previous work the experimental effort 
concentrated on the study of capture reactions using the 
activation technique and in a few cases on direct photo- 
disintegration studies with photon beams. The activa- 
tion technique for capture reactions on stable nuclei is 
limited to cases where the final reaction products are 
characterized by specific 7-radiation signatures. This 
is necessary for uniquely identifying the reaction prod- 
ucts and separating them from the large background 
activity originating from target impurities. This tech- 
nique for example can be applied for the measurement 
of 75 As(p,n) 75 Se, 106 Cd(a, 7 ) 110 Sn, or 196 Hg(a, 7 ) 200 Pb, 
as well as a larger set of (n,7) cross sections. 

If the characteristic activity is masked by background 
or if the reaction product lacks characteristic activity, al- 
ternative detection methods can be envisioned through 
AMS methods where the reaction products are chemi- 
cally separated from the activated target and analyzed 
through high resolution accelerator mass spectroscopy. 
This approach has been successfully appli ed in the study 
of th e 62 Ni(n,7) 63 Ni s-process reaction (|Nassar et al. I 
2005). Possible examples for applying this specific tech- 
nique are the reactions 72 Ge(p,7) 73 As, or for the higher 
mass range 156 Dy(a, 7) 160 Er or 190 Pt(a, 7) 194 Hg, which 
are characterized by such a low value that a measure- 
ment of the decay activity is difficult. 

Inverse photodisintegration measurements have 
also been developed as a powerful exp erimental tool 
l|Utsunomiva et al. I l200ll IVoet et al. I l200l(l . The 
increasing availability of high energy photon beams at 
the HIGS facility at TUNL (http://higs.tunl.duke.edu/) 
or at the ELBE facility of the Forschungszentrum 

Rossendorf 

( http : / / www . fz-rossendorf . de/pls/rois/ Cms?pNid= 144) 
opens new opportunities for targeting reactions such 
as 74 Se( 7 ,p) 73 As, 92 Mo( 7 ,p) 91 Nb, and 96 Ru( 7 ,p) 95 Tc, 
which directly affect the light p-nuclei abundances. 
Again, the reaction product can be detected by measur- 
ing the characteristic activity. If the latter is masked by 
background or falls below present detection limits alter- 
native analytical methods such as the AMS approach 



can be applied. 

Many of the reactions listed in Tables El and how- 
ever, are photodisintegration processes on radioactive nu- 
clei leading to radioactive nuclei; reactions with signifi- 
cant impact on flux and abundances are 110 Sn(7,p) 109 In, 
126 Ba(7,p) 125 Cs, or in the range N >82, for example 
152 Dy( 7 ,a) 148 Gd, 160 Yb( 7 , a) 156 Er. Such reactions can 
be approached by Coulomb dissociation techniques with 
radioactive beams. In the specific p-process cases listed 
here heavy ion radioactive beams have to be developed in 
the energy range of 2-12 MeV/amu depending on the as- 
sociated Q-values. Coulomb dissociation techniques are 
typically applied for studies with light radioactive par- 
ticles, which are produced through beam fragmentation 
reactions. Such cases are suited for facilities such as the 
NSCL at MSU, at RIKEN, and at GSI where heavy neu- 
tron deficient radioactive particles can be produced by 
spallation or heavy ion evaporation reactions at inverse 
kinematics. 

5. SUMMARY 

We have investigated the reaction flux patterns of the p 
process within the framework of a multi-mass zone type 
II supernova shock front model. Particular attention was 
paid to the influence of p-process reaction rates within 
the uncertainty limits of the theoretical Hauser- Feshbach 
model. The predicted p abundances are similar to the 
ones observed by previous work in this field. It could 
be shown that the endemic underproduction of the light 
p nuclei 92 ' 94 Mo and 96 ' 98 Ru are not due to the uncer- 
tainties in the thermonuclear reaction rates and must be 
traced back to other reasons. Possible neutron poison 
reactions in core helium burning and secondary neutron 
sources during core carbon burning could modify signif- 
icantly the weak s-process seed abundance distribution 
for the p-process. This was not considered in the present 
paper. It should be addressed in future studies to inves- 
tigate systematically the s-process related uncertainties 
for p-process simulations. 

Since (7,n) reactions control the overall feeding of the 
p-process flux, the corresponding reaction rates affect 
the final p-process abundance distribution over the en- 
tire mass range. This influence is particularly visible in 
the range above N =82. On the other hand we found 
that the impact of ( 7 ,p) and (7, a) reactions is limited 
to specific mass regions. Changes in the ( 7 ,p) reactions 
translate to direct modifications of the p-process abun- 
dances in the lower mass range with Z <50, but similar 
correlations were not observed in the higher mass range. 
In contrast, variations in the (7, a) reaction rates im- 
pact the abundance predictions in the higher mass range 
above ./V =82, while only small effects can be observed 
at lower masses. The overall impact of reaction rates 
seems not as dramatic as for other processes, the here 
discussed modifications of the rates translated in gen- 
eral to a change in the absolute abundances of a factor 
two to three. This will directly affect the overproduction 
factors displayed in Fig. Improved data are neces- 
sary since a reduction of the associated uncertainty will 
help to identify more clearly other uncertainty factors in 
the present simulation. This includes the uncertainties 
resulting from the initial seed distribution or uncertain- 
ties yielding from insufficient model descriptions of the 
p-process scenario. 
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Past experiments have confirmed that theoretically 
predicted p-process reaction rates agree within a factor of 
two with experimentally determined proton and neutron 
capture and the respective photodisintegration data. For 
the corresponding a rates considerably larger differences 
of up to a factor of ten were reported. More experimental 
work is necessary to test the validity of the a potentials 
used in the model calculations. It is also necessary to ex- 
pand the experimental work towards p-process reaction 
studies on neutron deficient nuclei to verify the reliabil- 
ity of the theoretical predictions for reactions involving 
radioactive nuclei. 
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TABLE 1 

MASS COORDINATE, MAXIMUM TEMPERATURE AND MAXIMUM DENSITY IN THE INVESTIGATED p-PROCESS layers. 



layer 
No. 


mass inside the shell maximum temperature 
(M ) (T 9 ) 


maximum density 
10 5 (g/cm 3 ) 


1 


1.9336 


3.45 


7.85 


2 


1.9658 


3.11 


6.64 


3 


2.0085 


2.96 


5.68 


4 


2.0508 


2.76 


4.82 


5 


2.1037 


2.60 


4.07 


6 


2.1564 


2.44 


3.56 


7 


2.2090 


2.32 


3.16 


8 


2.2614 


2.21 


2.80 


9 


2.3136 


2.12 


2.54 


10 


2.3655 


2.04 


2.30 


11 


2.4171 


1.97 


2.11 


12 


2.4684 


1.91 


1.94 


13 


2.5249 


1.84 


1.75 


14 


2.5825 


1.79 


1.68 
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TABLE 2 

Selected (7,p) or (n,p) reactions that together with their respective inverse reactions were fou nd t o exhibit the 
strongest influence on the final p abundances. their impact is illustrated in the top panel of flg. ^] particularly 

important rates are marked with *) . 



126 Ba( 7 ,p) 125 Cs*> 92 Mo( 7 ,p) 91 Nb*-> 
110 Sn(7,p) 109 In*> 86 Rb(n,p) 86 Kr*) 



106 Cd( 7 ,p) 105 Ag 
104 Cd( 7 ,p) 103 Ag 
100 Pd( 7 ,p) 99 Rh 
96 Ru( 7 ,p) 95 Tc*> 



85 Sr(n,p) 85 Rb*) 
84 Sr(7,p) 83 Rb*) 
78 Kr( 7 ,p) 77 Br*) 
77 Se(n,p) 77 As 



75 So(n,p) 75 As*) 

74 Se( 7 ,p) 73 As*) 

76 As(n,p) 76 Ge*> 

75 As( 7 ,p) 74 Ge 

73 As( 7 ,p) 72 Ge*) 

71 Ge(n,p) 71 Ga 
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TABLE 3 

Selected (7, a) reaction chains, which were found to exhibit the strongest influence on the final p abundances. Their 

IMPACT IS ILLUSTRATED IN THE BOTTOM PANEL OF FlG. ^] PARTICULARLY IMPORTANT RATES ARE MARKED WITH 



74 Se( 7 ,a) 70 Ge*) 

96 Ru(7,a) 92 Mo*> 

106 Cd( 7 ,a) 102 Pd 

110 Sn( 7 ,a) 106 Cd*) 

120 Te(7,a) 116 Sn*> 

122 Xe( 7 ,a) 118 Te 

128 Ba( 7 ,a) 124 Xe*) 



196 Pb( 7 ,a) 192 Hg*> 

191 Hg(27,2a) 183 Os 

190 Hg(37,3a) 178 W*) 

189 Hg(3 7 ,3a) 177 W 

188 Hg(57,5Q) 168 Yb*) 

183 Pt(37,3a) 171 Hf 

178 Os(4 7 ,4a) 162 Er*> 

177 Os(3 7 ,3a) 165 Yb 

176 Os(57,5a) 156 Dy*) 

167 Hf( 7 ,a) 163 Yb 

166 Hf(5 7 ,5a) 146 Sm*) 

156 Er(37,3a) 144 Sm*> 



Sensit 

io 2 n 




TIME (s) TIME (s) 

Fig. 2. — Temperature and density profiles of the supernova shock front traversing the Ne/O layer of the pre-supernova star. 
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Fig. 3. — Adopted seed abundance distribution for the p-proccss simulations 



100 125 150 

MASS NUMBER 



175 



200 



Sensitivity of p F 



>va Model 17 



A 
V 




120 140 
MASS NUMBER 



200 



Fig. 4. — Averaged normalized overproduction factor for the proton rich p nuclei from network calculation with standard reaction rates 
(see text). 
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Fig. 5. — Integrated reaction flux during the first second of a type II supernova explosion in the Ne/O layer with a maximal temperature 
of T9=2.96. The reaction flux is indicated by the line thickness (thick solid line >10~ 10 , solid line >10 -11 , dashed line >f0~ 12 ). 



Sensitivity of p Process Nucleosynthesis to Nuclear Reaction Rates in a 25 Solar Mass Supernova Model 19 




75 77 79 81 83 85 87 89 91 93 95 97 



Cd (48) | 





















£ 

lb (65) 


yj.66) 


























1 1 1 1 




































Xe 
e(52) 


54) 


Ba (5f 
Cs (55) 


5) 


L 


Gd (64) 










































99 


Eu 


63) 












































Sm (62) 




















































Pm (6 


















































Nd (60) 
























































Pr 

Ce (5 
a (57) 


59) 














































































































F 


















































yo 










































































































































u 




































I (53 






















































oy y i yo 


























































Of 


Sb 


51) 
















































J) 






















































85 
































































rhr 






























81 82 





53 55 57 59 61 63 65 67 69 71 73 75 77 79 




31 33 35 37 39 41 43 45 47 49 



Fig. 6. — Q-values for (7,11) reactions based on experimental mass data lAudi et al. ll200l) . 
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Fig. 9. — The integrated reaction flux during the first second of a type II supernova explosion in the Nc/O layer with a maximal 
temperature of T9=2.44. The reaction flux is indicated by the line thickness (thick solid line >10 -11 , solid line >10~ 12 , dashed line 

>io- 13 ). 
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smaller and larger rates, respectively. 
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Fig. 11. — Overproduction factors for selected p nuclei ( 96 Ru, 98 Ru, 112 Sn, 138 La, 144 Sm, 162 Er and 184 Os) as a function of temperature 
indicating the production efficiency in the different p-process zones. 
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Fig. 12. — Same as bottom panel of Fig. 1101 but with the (7,0) rates modified by factors of 0.1 (solid line) and 10 (dashed line). The 
results for three times smaller and larger rates are indicated by squares and crosses, respectively. 
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Fig. 14. — Ratio of the overproduction factors for the p abundances in the mass region 92< A < 98 obtained when reducing the (7,p) 
destruction rates of 92 Mo, 94 Mo 96 Ru, and 9s Ru by a factor of ten to the overproduction factors obtained with unchanged rates. 
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Fig. 15. — Influence of the selected reactions listed in Tables 01 and 131 compared to the results obtained by modification of all accepted 
HF-rates as shown in Fig. 1101 ('symbols'). The top panel refers to the proton reactions of Table 01 and the bottom panel to the a reactions 
of Table |3 The effect of chan ges in the selected reactions is indicated by solid and dashed lines for reducing and increasing the rates by a 
factor of three, respectively. It is obvious that the selected rates are almost completely accounting for the response of the p-process network 
to reaction rate uncertainties. 
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Fig. 16. — Cross section data for two sets of 75 As(p,n) me asurements ((squares ) IKaila^et al. 1119791) . (circles) IMushtaaet et al. 11198a) 1 
in comparison with Hauser Feshbach predictions (solid line) IHauscher fc Tnielemamil20'0l|) . Also shown is the experimental cross section 
of 85 Rb(p,n) ((crosses) IKastleier et al. 11200^1 1 in comparison with Hauser Feshbach predictions (dashed line) (Rauschcr & Thiele mann I 
2001). 



